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A preliminary weather model is presented for optical propagation through the atmo- 
sphere . It can be used to compute the attenuation loss due to the atmosphere for desired 
link availability statistics. The quantitative results that can be obtained from this model 
provide good estimates for the atmospheric link budget necessary for the design of an 
optical communication system. The result is extended to provide for the computation of 
joint attenuation probability for n sites with uncorrelated weather patterns. 


I. Introduction 

During the past few decades, the pioneering work of 
Tatarski [1], Fried [2,3], Hufnagel and Stanley [4], and 
Ishimaru [5] , has set the stage for the development of a con- 
sistent theory to evaluate the performance of optical commu- 
nications and imaging through the atmosphere. However, there 
is a conspicuous lack of weather models that may be used to 
obtain even a first-order estimate of statistics for the avail- 
ability of an optical communications link under ambient 
weather conditions. 


The light energy of a beam dissipates as it travels through 
the atmosphere, due to scattering and absorption. If it can be 
assumed that (1) attenuation loss is independent of radiation 
intensity and (2) the absorbing and scattering events occur 
independently of each other, the atmospheric attenuation can 
be expressed by the Bouguer-Lambert law: 


Iiy) = I Q {v) exp 



y(v,z)dz 


( 1 ) 


where I(v) is the observed irradiance at optical frequency v, 
/ 0 (v) is the irradiance that would have been observed if the 
beam were propagating in a vacuum at a distance Z from the 
source, and y(v, z ) is the total extinction coefficient due to 
scattering and absorption by the atmospheric constituents at 
position z. The magnitude of the argument of the exponential 
in the above equation is known as the optical depth or optical 
thickness, r, of the medium, i.e.. 


T 



y(v, z)dz 


( 2 ) 


Experiments [6] with artificial fog and smoke and diluted 
solutions of milk show that the Bouguer-Lambert law holds 
well for optical thickness r < 12. The optical depth can be 
simply related to the attenuation loss, L, in dB by the fol- 
lowing approximation: 


L »» 4.34r 


( 3 ) 
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The distribution of gas molecules, clouds, fog, haze, aero- 
sols, and other particulates in the transmission path influence 
the computation of attenuation loss. These phenomena at best 
can be interpreted on a statistical basis. The system designer 
may have to be content with a probability for which the atten- 
uation of the atmosphere due to ambient weather is less than 
some critical value Z dB. Hence we define the probability, 
w^Z), of the random weather variable, W , for a single site as 

Wl (L) = P W (KL) (4) 

where / is the loss variable. wq(Z) is the fraction of weather 
conditions at any given site for which the attenuated direct 
beam irradiance /(v) > / Q (v) exp [-0.23 Z] . The weather 
model developed below can be used to compute this proba- 
bility. The result is also extended to obtain the joint weather 
probability, h^(Z), that at least one of the sites has extinction 
loss / < Z forn independent sites receiving simultaneously. 


rapidly with increasing concentration of aerosols, fog, 
haze, and clouds. Optical attenuations L > 1000 dB are 
possible, where the validity of the Bouguer-Lambert 
law is questionable. Table 1 shows typical ranges of 
attenuation values for various types of clouds [8] . 
Since it is unlikely that an optical communication 
system will be designed for Z > 30 dB, the Bouguer- 
Lambert law serves as a good approximation. As shown 
in Fig. 2, the probability of opaque cloud cover over 
the southwestern U.S. is 20% [7] . Let us take a pessi- 
mistic view and assume that attenuation loss due to all 
opaque clouds is 100 dB or higher, and use this assump- 
tion to estimate one of the model parameters. 

(4) It is further assumed that the attenuation of the beam 
has an exponential distribution for L >L 0 . The hypo- 
thesis is supported by visual observations, and also by 
the experimental data collected at Goldstone in Cali- 
fornia for atmospheric propagation at 8 to 10 Ghz. 1 


II. Weather Model 

The following observations are made in order to arrive at 
the preliminary weather model: 

(1) It is assumed that for some fraction of the time p, 
where 0 < p < 1, clear weather conditions hold. This 
number can be determined approximately from exist- 
ing data on cloud cover and visibility for potential 
sites. An analysis of two years of GOES satellite data 
by Wylie and Menzel [7] shows that the probability 
of having clear weather in the southwestern U.S. is 
over 60% (Fig. 1). 

(2) The attenuation under clear weather conditions is due 
to scattering and absorption by air molecules and 
sparse particulate matter in the upper atmosphere. .The 
attenuation due to molecular absorption is approxi- 
mately 0.5 dB [6] . The attenuation due to molecular 
scattering is of the same order. Given that the atmo- 
sphere is never totally free of aerosols and thin cirrus 
clouds, on an average clear day the attenuation loss 
would be in the range of 1 to 3 dB. For simplicity, 
let us assume that the minimum attenuation loss due 
to the atmosphere is 3 dB. Hence, the model defines 
clear-air atmosphere, which occurs with probability 
p, as having an attenuation lossZ 0 = 3 dB. Note that 
this attenuation loss refers to near-zenith propagation 
through the entire atmosphere. For any zenith angle 0 , 
the path attenuation can be written as Z Q sec 6. How- 
ever, the remainder of the discussion assumes near- 
zenith propagation paths, i.e., 0=0. 

(3) The probability of not having clear weather isq = 1 -p. 
The attenuation loss of the atmosphere increases 


With the foregoing assumptions in mind, it is now possible 
to postulate a weather model. The weather cumulative distri- 
bution function (CDF) for a single site can be modeled as 

w'jCO = 1 -?exp [-0.23 *(Z. - i 0 )] {L>L Q ) 


( 5 ) 

where w^Z) is defined in Eq. (3) above, q is the probability 
fraction when the weather is not clear, b is a site-dependent 
parameter to model the slope of the CDF curve, and Z 0 is 
the minimum attenuation loss due to the atmosphere. Note 
that w x (L = L 0 ) = p . The value of b may vary with geo- 
graphical location and altitude, and can be inferred from 
observed visibility and extinction loss data for potential 
receiving sites. For the southwestern U.S. region, q = 0.4, 
and from the third item above, the attenuation loss is 100 dB 
when w^L = 100) = 0.8. Using these values in Eq. (5), we find 
that b = 0.03. This estimate of the parameter b will be used in 
numerical examples later. 


Equation (5) can be recast in a more familiar form for an 
optical communication system designer to give the attenua- 
tion loss in terms of the weather probability uq(Z), i.e., 


L = 


L o . 


Z rt + 


1 


o 0.23 b 


In 


1-VZ) 


for vv 1 (Z) <p 
for (Z)> p 


( 6 ) 


1 JPL Internal Document 810-5, Rev. D, 1988. 
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Figure 3 is a plot of attenuation loss as a function of the 
weather probability for the southwestern region of the U.S., 
for which q = 0.4, b =0.03, and L 0 = 3 dB. These estimates, 
needless to say, are quite important for the system designer to 
determine the link budget for loss due to the atmosphere for 
Earth-space optical communication paths. 

It is also possible to extend the result to n sites receiving 
simultaneously. From Eq. (4), it is easy to see that the comple- 
ment of the probability w^Z,) is given by 

P w (l>L) = w'd) = qexp[-0.23b(L-L 0 )] 

(L>L 0 ) (7) 

The joint probability of n sites, vv£(Z,), that the attenuation 
loss /• > L for all i can be written as 

<a) = „„(/,> l,i 2 > l,..., / n >D 

( 8 ) 

where the subscripts 1 , 2, . . . , n label the receiving sites. When 
the weather conditions for all sites are independent and iden- 
tically distributed (IID), we have 

w„ c (Z) = P Wi (l t > L) P W2 (l 2 P w Jl n > L ) 

= [P w U>l)] n (9) 

Using Eq. (9), the joint probability that at least one site has 
attenuation KL is found to be 


w n (L) = !-<(/,) = 1 - [q exp [-0.23 b{L-L Q )]] n 

(L>L 0 ) (10) 

For a single site with p = 0.6, L Q = 3 dB, and b = 0.03, the 
probability that the attenuation/, < 3 dB is = 3) = 0.60. 
If three such IID sites are chosen, we have w 3 (Z = 3) = 0.94. 
In other words, if a system is designed to absorb an extinction 
loss of 3 dB, a three-site receiving network will be functional 
94% of the time. Figure 4 plots the fraction of the total period 
under ambient weather conditions when the attenuation is 
< 3 dB as a function of the number of sites. Table 2 gives the 
expected dB loss for a desired link availability percentage for 
up to four joint receiving sites. It is, however, not very clear 
how the independence of weather patterns at various sites can 
be insured. It is known that the scale size of weather patterns 
is on the order of a few hundred kilometers, and this measure 
may be used to find sites with uncorrelated weather. Joint ob- 
servation of weather parameters for the probable sites will be 
necessary to make a more accurate determination. 

III. Conclusion 

The virtue of the weather model presented here lies in its 
simplicity. It may be applied with ease to obtain a first-order 
magnitude of probabilities w n (L) for an optical system that 
must operate in the atmosphere. The computation of these 
probabilities, for example, will provide a good statistical esti- 
mate of the attenuation loss to an optical communication 
system designer for link availability. 

The model does not consider frequency dependence, since 
it has been studied thoroughly and well documented in LOW- 
TRAN computer code developed at the Air Force Geophysics 
Laboratory [9] . The model also disregards seasonal variations, 
which can be incorporated later when adequate data bases 
have been developed. 
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Table 1. Typical range of values for cloud attenuation 
(adapted from [8]) 


Cloud 

type 

Cloud base 
height, 
km 

Cloud 

thickness, 

km 

Total 

vertical 

attenuation, 

dB 

Stratus 

0.1 - 0.7 

0.2 - 0.8 

26 - 350 

Stratocumulus 

0.6 - 1.5 

0.2 - 0.8 

3 - 100 

Nimbostratus 

0.1 - 1.0 

2 - 3 

260 - 1300 

Altostratus/ 

cumulus 

2-6 

0.2-2 

9 - 260 

Cumulus 

0.5 - 1.0 

0.5 -5 

22 - 870 

Cumulonimbus 

0.5 - 1.0 

2-12 

260 - 5200 

Cirriform (ice) 

6-10 

1.0 -2.5 

1 - 15 

Fog 

0 

0-0.15 

0-13 


Table 2. Attenuation loss as a function of desired link 
availability for n sites receiving jointly 

Percentage 

weather 


Attenuation loss, L y dB 


n = 1 

n = 2 

n = 3 

n = 4 

60.0 

3.00 

3.00 

3.00 

3.00 

62.0 

10.34 

3.00 

3.00 

3.00 

64.0 

18.07 

3.00 

3.00 

3.00 

66.0 

26.25 

3.00 

3.00 

3.00 

68.0 

34.92 

3.00 

3.00 

3.00 

70.0 

44.16 

3.00 

3.00 

3.00 

72.0 

54.03 

3.00 

3.00 

3.00 

74.0 

64.63 

3.00 

3.00 

3.00 

76.0 

76.08 

3.00 

3.00 

3.00 

78.0 

88.53 

3.00 

3.00 

3.00 

80.0 

102.16 

3.00 

3.00 

3.00 

82.0 

117.23 

3.00 

3.00 

3.00 

84.0 

134.09 

3.00 

3.00 

3.00 

86.0 

153.19 

12.55 

3.00 

3.00 

88.0 

175.24 

23.58 

3.00 

3.00 

90.0 

201.32 

36.62 

3.00 

3.00 

92.0 

233.25 

52.58 

3.00 

3.00 

94.0 

274.40 

73.16 

6.08 

3.00 

96.0 

332.41 

102.16 

25.41 

3.00 

98.0 

431.57 

151.74 

58.47 

11.83 

99.0 

530.72 

201.32 

91.52 

36.62 
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PROBABILITY OF CLEAR SKY, PERCENT 

Fig. 1 . Contour diagram obtained from 2 years of GOES satellite data showing the probability 
of clear sky over the United States [7]. 



PROBABILITY OF OPAQUE CLOUD COVER, PERCENT 

Fig. 2. Contour diagram obtained from 2 years of GOES satellite data showing the probability of 
opaque cloud cover over the United States [7]. 
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